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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 28^9 

CORRECTIONS FOR LIFT, DRAG, AND MOMENT OF AN AIRFOIL 
IN A SUPERSONIC TUNNEL HAVING A GIVEN 
STATIC PRTSSURE GRADIENT 
By H. F. Ludloff and M. B. Friedman 

SUMMARY 


The corrections for lift, drag, and moment of a two-dimensional air- 
foil have "been analyzed, on: the assumption that the airfoil is tested in 
the working section of a supersonic tunnel in which the pressure field, 
instead of being uniform, is characterized by gradients in the axial and 
transverse directions. 

The pressure gradients of the tunnel as well as the effect of the 
airfoil to be tested are regarded as perturbations of the original rec- 
tilinear flow field of given Mach number. Therefore the velocity poten- 
tial of the flow, the nonlinear differential equation of motion, and the 
boundary conditions are expanded into double series in powers of two 
parameters, one characterizing the airfoil thickness e and the other, 
the inhomogeneity of the field b. In this way the nonlinear problem 
is split into a system of linear boundary- value problems, corresponding 
to the different powers of b and e. 

In each of the resulting problems there appears, besides the dif- 
ferential equation and boundary condition, an additional condition to 
be stipulated on the characteristics passing through the leading edge. 
Particular attention has been paid to the correct formulation of this 
"characteristic condition." 

The solution procedure is carried out up to orders b^, e^, and 

eb in the velocity potential. This means that, for example, the drag 
is computed up to orders b^e, , and e^b. The physical meaning of 
the results is discussed. The drag term in eb represents the 
"horizontal buoyancy" of the airfoil, the term proportional to be^ is 
a consequence of the interaction of the airfoil field and the Inhomo- 
geneous pressure field, and the term in b^e may be considered as a 
"second-order buoyancy." The meaning of the various lift and moment 
terms may be interpreted similarly. The resulting expressions have been 
derived for arbitrary given pressure gradients and general profile form. 

All solutions are obtained in closed, analytic form ready for imme- 
diate evaluation. Representative examples with graphs are included. 
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INTRODUCTION 


Reports from supersonic tunnels indicate that sometimes there exist 
in the working section undesired static pressure gradients due to design 
or construction flaws, which are difficult to eliminate. As a consequence, 
force measurements made upon models have to he corrected to agree with the 
forces to he expected in rectilinear flow. 

For incompressible flow the correction to he applied to the drag, 
in the case of a longitudinal pressure gradient is well-known (refer- 
ence l). Because of the occurrence of a "horizontal buoyancy, " the cor- 
rection turns out to he equal to the product of the pressure gradient 
and the sum of volume and "apparent volume" of the test body; the appar- 
ent volume, like the real volume, depends only on the geometrical shape. 

The present analysis is an attempt to solve the analogous problem 
in supersonic flow under the most general conditions, which require cor- 
rections to he applied also to the lift and moment. A too -dimensional 
pressure gradient in both the longitudinal and transverse directions is 
considered, which is equivalent to a stream-angle variation along the 
tunnel axis. In the first approximation the computations must he 
expected to yield a superposition of well-known results. Hence, it is 
important to carry the calculations up to the second approximation, 
which includes a term characterizing the interaction between the airfoil 
field and that of the given pressure gradient. 

This work was conducted at New York University under the sponsorship 
and with the financial assistance of the National Advisory Committee for 
Aeronautics . 


SYMBOLS 


h parameter characterizing inhomogeneity of field 

c local sound speed 

c Q sound speed at Mq 

k(x) dimensionless airfoil profile function measured 

above flight direction 

k(x) profile function measured above chord 

Z * chord length 
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Po 

q 

p = \[ Mo 2 - 1 
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e 

e 

e 

p 

p o 

9 

9 s 

9 ' 

t(x) 


local Mach number 

Mach number of uni form field 

x- and y-components of velocity, on axes, due to 
inhomogeneity of flow field (dimensionless) 

o 

local pressure 
pressure at Mq 
local velocity 

dimensionless coordinates, axial and vertical 


adiabatic exponent 

airfoil thickness parameter i 

local angle of attack 

angle of attack of chord 

characteristic coordinates 

local density 

density at Mq 

dimensionless velocity potential 
velocity potential of basic flow field 
velocity potential due to airfoil disturbance 
ordinate of leading-edge characteristic 


STATEMENT OF PROBLEM 
Physical Formulation 


Suppose that in the test section of a supersonic tunnel the velocity 
field, instead of being uniform, is represented by the velocity distribu- 
tion along the axis y = 0, as follows: 
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At y = 0: 


Mx = Mq + bm(x) 
My = bn(x) 


( 1 ) 


Here M x = q x /c Q and My = qyJc Q} where q x and qy are the x- and 
y- components of the local velocity of the flow field, c Q is the sound 
velocity in a region where the Mach number is Mq, x and y are 
dimensionless coordinates measured in units of the length l of the test 
body, and b is a dimensionless parameter characterizing the magnitude 
of the inhomogeneity of the flow field and the pressure gradient. Under 
actual conditions b is a small quantity such that 


bm(x) 

bn(x) 


« Mo 


( 2 ) 


where m(x) and n(x) are given arbitrary functions along the axis, 
which may be determined by measurement. With conditions on the axis 
fixed the flow field in the whole xy-plane is automatically prescribed 
by the equations of supersonic flow. Terms of the order b 2 may also 
be given arbitrarily on the axis but are assumed to be zero in the 
present case. This does not imply, however, that terms of the order b 2 
do not appear in the field away from the axis. 

In this velocity field a thin airfoil of given profile form k(x) 
and thickness ratio e is inserted, which may be regarded as another 
small perturbation of the rectilinear flow field Mq. Two kinds of dis- 
turbances then exist, the one produced by the velocity or pressure gra- 
dient measured by b and another one produced by the airfoil and 
measured by e. Therefore, it appears appropriate to expand the veloc- 
ity potential cp, as well 'as the potential equation of supersonic flow, 
and the boundary condition on the airfoil surface into a power series in 
powers of both parameters b and e . , 

In this way, the nonlinear problem is split into a system of linear 
boundary-value problems corresponding to the different powers of b 
and e. In each of the resulting problems there appears, besides the 
differential equation and boundary condition, an additional condition to 
be stipulated on the characteristics passing through the leading edge: 

It is to be required that the disturbance produced by the airfoil vanish 
along these curved lines in order to make the solution unique. It may 
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"be noted that the characteristics of the complete differential equation 
are curved, while those of the resulting linear equations are straight. 
It is shown later how the condition on the characteristics can be for- 
mulated for each linear boundary- value problem separately. 


Analytic Formulation 

Because of the relatively small thickness of the airfoil, the shocks 
at the leading edge will be weak and the entropy changes in the flow 
negligible. Hence, the complete flow field will be described by a veloc- 
ity potential which obeys the well-known, nonlinear differential equation. 
Introducing a nondimens ional quantity cp, equal to- the potential divided 
by Zc Q , this equation may be written: 



<Px 2 ) c Px X - ^9x < Py t Pxy + i~2 " = 0 


(3) 


where the subscripts x and y indicate differentiation. 

The local sound speed c corresponding to the local velocity q 
is related to the sound speed c Q by means of the energy equation 


where 



q 


2 





(4) 


Substituting equation (4), equation (3) can be written: 


E + - •h 2 ) 

■ <fx2 ) 


4^ %T\ lixx - + 


T 1 ^3 V ■ 0 


2 


(5) 
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Row, 9 may be expanded as a series in powers of the perturbation 
parameters b and e, an follows: 

9 = cp 00 + bcfp-*- + b ^ 02 + bcp^ - ^ + b 2 cp^^ + ecp ^ 2 + ebcp^ + £2^4 ( 6 ) 


The potential coefficients 

cp 8 = cp 00 + b9 01 + b 2 9 02 (7) 

represent the potential of the basic inhomogeneous velocity field before 
an airfoil was placed in it. 

The terms 


9' = b9 1< - > + b 2 9 1 - 1 - + e 2 9- 1 - 2 + ebqr 1 -^ + g^cpl^ - (8) 


represent the perturbation potential produced by the airfoil in the 
basic field. 

It is to be noted that the perturbation potential contain^ the 
terms bqp 10 +' b'vp 11 which do not vanish for an airfoil of zero thick- 
ness ( e = 0) . 

A disturbance for which e = 0 may be thought of as a very thin 
flat plate placed at a given angle into the basic flow field. In order 
for the inhomogeneous flow to follow the surface of the plate, vertical 
disturbance velocities must be produced of such magnitude (up or down- 
ward) as to cancel the vertical velocity components of the inhomogeneous 
flow field, which are proportional to b or b 2 . 

The potential term proportional to eb is due to the interaction 
between the basic field and that produced by the airfoil. The e and 
e 2 terms represent the first- and second-order potentials of an airfoil 
in a homogeneous field. 

Finally, the analytic procedure to be followed in this investigation, 
will be outlined briefly: 

(l) Differential equation: Substituting series (6) for 9 into the 

nonlinear equation (equation (5)) and splitting in powers of b and ■£, 
a set of homogeneous or nonhomogeneous equations of the wave type is 
obtained. 
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Those equations which are proportional to b, e, and are 

homogeneous : 

(Mo 2 " ~ <Pyy lk = 0 (9) 

O 

The equations which are proportional to b and eb are of the 
nonhombgeneous type : 

K 2 - <Pyy ik = *±k{v lm ) (10) 

where the functions can be shown to be known functions of lower- 

order cp im which have been computed before. 

(2) Boundary condition: Insert series (6) for cp into the usual 

boundary condition stating that the flow is everywhere tangential to the 
surface of the body: 

l~] = ek* (x) (11) 

yPx /surface 

and expand equation (ll) about y = 0 in powers of e. This imposes 
on each of the potential coefficients cpl^ a condition of the form 

(V% 0 ■ °^ lm ) <12) 


where, again, the functions G-y^- can be shown to be known functions of 
lower-order cp im . 

(3) "Characteristic condition": As mentioned before, the poten- 

tial cp 1 describing the various disturbances produced by the airfoil 
is required to vanish on the leading-edge characteristics which may be 
written: 


y = i|r(x,b) = ->Jr 0 (x) + bi^x) + . . . (13) 

since the form of the characteristics is influenced by the inhomogeneity b. 
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It can be shown (see section "Characteristic Condition") that the 
"characteristic condition" for cp* is equivalent to the following 

11c 

relations to be imposed on the individual <p : 

2 

For the orders b, e, and e 


<^(*,*0) = ° 


( 14 ) 


For the orders 


and be 


cp^x,^) = i|r 1 (x)cp 2m (x,i|ro) (15) 


where , again, the coefficients q)^ m are known lower-order potentials 
and Tjf Q (x) and t|/j(x) may be computed by means of the theory of 

characteristics . 

At the end of this procedure the solutions <p*^ will have been 
derived in closed, analytic form, and the pressures at any point can be 
computed. Since the objective is to compute the corrections for forces 
and moments, it is sufficient to determine the velocities and pressures 
on the airfoil surface only, thus simplifying the considerations. 


BASIC VELOCITY FIELD 
Differential Equation 


The analysis can be clarified by considering first the basic inhomo- 
geneous field alone. Then the series expansion for cp reduces to: 


qp® = qP® + bqp®l + b^cp 02 


(16) 
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Substituting equation (l) into equation (5) yields: 

|l - Mo 2 - b(y + l)Mo<Px 01 + • * ^]( b %cx 01 + ‘b 2 ^ 02 +•••)- 
2(Mq + . . ^(bqy 01 + . . .)(b V 01 + . . . j + 

[l - b( 7 - DM^ 01 + . . (bcPyy 01 + ^ 2 <tVy° 2 +•••)= 0 (17) 

In deriving equation (17) cp 00 = MqX has been substituted in equa- 
tion (5), since it is obviously a solution; furthermore, all those terms 
have been kept which could possibly contribute to the differential equa- 
tion in b and b 2 . Splitting up in powers of b, the following equa- 
tions result: 

In b: 


K 2 - 1 )i> xx 01 



= o 


p 

and in b : 


(Wo 2 - i) T <„ 02 - Tyy 02 - -MojTr + D^ 0 ^ 01 + (7 - 1 H> x 0l 9 yy 01 + 

w*] 


Boundary Condition 

The velocity distribution (equation (l)), prescribed on the axis 
y = 0, may be considered as a boundary condition for the potential of 
the basic flow. 
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and 


Substituting equation (l6) into equation (l), it follows that 

O' 

Mo + b(p x 01 (x_,0) + b 2 qp x 02 (x,0) = Mq + bm(x) 
bcp y 01 (x,0) + b 2 cp y 02 (x,0) = bn(x) 


(18) 


Ordering terms in powers of b, a boundary problem in b and one in b 2 
are obtained. 


Boundary- Value Problem in Order b 

The boundary- value problem in order b consists of the differential 
equation 


p 2 qp 01 - cp 01 = 0 
H T xx Y yy 


where p 2 = Mq 2 - 1, and the boundary conditions, at y = 0, 


(19a) 


cp x 01 (x,0) = m(x) 


qy^x^O) = n(x) 


(19b) 


The solution of differential equation (l9a) is 

Cp 01 (x,y) = f(x + py) + g(x - py) 


( 20 ) 


where f and g are arbitrary functions to be determined by boundary 
condit ions ( 19 b ) . 

Differentiating equation (20) with regard to x and y and then 
putting y = 0 yield: 


01 / 


( 21 ) 


<Px (x,0) = f» (x) + g’ (x) = m(x) 

cp y 01 (x,0) = p[f' (x) - g’ (xj| = n(x) 
where f 1 and g 1 indicate derivatives with regard t'o the whole argument . 


i 



KACA TN 2849 


11 


From equations (21) there follow: 



( 22 ) 


Since it can be seen from equation (20) that f and g or their 
derivatives are constant along their characteristics x ± py, it may be 
concluded from equations (22) that 


and 


f* (x + Py) = | 

m(x + py) + i n(x + py)J 

g* (x - py) = | 

m(x - py) - i n(x - py)J 


(23) 


Differentiating equation (20) and substituting equations (23), the 
velocity components of order b 1 for the entire field follow: 


9 x 01 (x,y) = |[m(x + Py) + | n(x + py) + m(x - p y ) _ i n (x 

<Py 01 ( x ;y) = |[™( x + Py) + | n(x + py) - m(x - py) + | n(x 

From equations (24) the expressions for cp^ 01 and cp xy 01 
in the differential equations for cp^ - ^ - and cp^*^ may be readily derived* 

Since, later on, the values of the velocity and pressure field are 
required only on the surface of the body, that is, near the axis, expres- 
sions (24) may be expanded about the axis: 

9 x 01 [x,ek(xi] = m(x) + ek(x)n’ (x) ... 1 



<Py 01 [x,ek(x)] = n(x) + ek(x)p 2 m’ (x) . . . 


( 25 ) 
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Boundary- Value Problem in Order b 2 


In order b^ 
ential equation 


tbe boundary- value problem consists of tbe differ- 


P 2 cp. 


02 


xx 


- 9yy 02 = + l)^ 0 ^ 01 + (7 - Dcp^CPyy 01 + 



( 26 a) 


and the boundary conditions 


OP 

<Px (x,0) = 0 

<fy 02 (x,o) = o 


( 26 b) 


In order to determine the velocity and pressure field on the airfoil sur- 
face, no explicit solution of equations (26) is necessary. This can be 
seen by the following reasoning. 


For y = ek(x) the derivatives of cpy02( Xj y) ma y be obtained by 
expanding about the axis: 


qj^^ekU)] = qj/^x^O) + ek(x)cp xy 02 (x,0) + 
q)y 02 [x,ek(x)] = qy 02 (x,0) + ek(x)cp yy 02 (x,0) + 


^ (27) 


Since according to equation (l6) cp^ 2 is of order b 2 , the second- and 
higher-order terms in equations ( 27 ) are of order b2g, b2g2 # . . 
and may thus be neglected so that in order b : 


<Px° 2 Jx,ek(xTl S 0 


q> y 02 [x,ek(x)] ~ 0 




(28) 
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Therefore, the total velocity components of the basic field on the 
airfoil surface are: 


cPx® = Mq + bm(x-) + 6bk(x)n* (x) + b 2 0 . . 

I 

<Py B = bn(x) + ebk(x)p 2 m t (x) + b 2 0 . . . 


( 29 ) 


It is to be noted that the second derivatives cp^® 2 and 9yy^ 

appear in the differential equation for <pll (see equation (32c)) and 
hence it may seem that the boundary- value problem (equations ( 26 )) has 

to be solved after all. However, the combination Tyy 02 - p 2 cp xx 02 , 
which is required in equation (32c), may be expressed in terms of the 
known lower-order potential <p°l, as differential equation ( 26 a) indi- 
cates. Carrying out the proper differentiations and substitutions 
yields: 



where £ = x + f3y and q = x — py are the characteristic coordinates 
of the problem. 



14 


NACA TN 2849 


VELOCITY FIELD OF AIRFOIL 
Differential Equation 


The next objective is to determine the perturbation potential pro- 
duced by an airfoil which is placed into the (inhomogeneous) basic field 
at an arbitrary small angle of attack, the leading edge coinciding with 
the coordinate origin. Inserting series (6) for cp into equation (5), 
the quantities cp 01 and qp^ 2 may now be considered as known and 

... cp 11 *, as looked for. Keeping only those terms which may pos- 
sibly contribute to a differential equation of the order b, e, b 2 , 
be, or e 2 , (equation 5) yields: 

|l - Mo 2 - 2L±J:|^2M 0 ((p x 01 + cp/ 0 ) + e2M Q ( ?x 12 . . 7|T- x 

17/ 01 10\ 12 ,2/ 02 ll\ 13 2 l4 “1 

[H^cx + <Pxx ) + e< Pxx + b (<Pxx + 9xx j + ebcp xx + e cp JfX . . .J - 


2 M o^(c P y 01 + cp/ 0 ) + ec/ 2 . . 7| [b(<p xy 01 + cp*/ 0 ) + ecp*/ 2 . . ^ + 



eqp. 


12 


yy 


W 2 


^yy 


ii 


+ eb 


™ry 


13 


+ = 


(3D 


Ordering terms in powers of b and e, equation (31) will split into 
the following set of linear equations: 

In b: 


□ 2 m 10 
P 'Pxx 


- cp. 


yy 


10-0 


(32a) 
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In e : 


In b 2 : 


V%oc 12 ~ ^ = 0 


(32b) 


P^x 11 - 'Pyy 11 = -%P + 1) + (7 - l)p!J(q> xx 01 + ^cx 10 )^ 01 + <Px 10 ) 


’(cpy 01 + c Py 10 )( t P X y 01 + ®_ 10 } U - 02 ° 2 ~ 02 


+ V AV + V ) t + Vyy - t v™ (32c) 


In be : 


P^xx 13 - V 13 = _M ° 


[(7 + 1) + (7 - l)pf] <P X 12 (<P ; 


01 + ^ 10 ) + 


XX T xx 


■fttx 12 ^ 01 + %t 10 )] + 2 j^ 12 ( , Pxy 01 + 'Pxy 10 ) + 

v 12 K 01 + V 0 )]} (32a) 


In e 2 : 


P 2( Pxx l4 - ( Pyy 1 ^ = -M 0 [7 + l)<P x 12 <Pxx 12 + (? ~ D'Px^'Pyy 12 


2cp -l- 2 © ^ 2 


(32e) 


Whenever the differential equation is inhomogeneous, the right-hand side 
turns out to be a function of the lower-order solutions which have been 
computed in an earlier step. 


Boundary Condition 

The usual boundary condition that the velocity of the air is every- 
where tangential to the airfoil surface may be written as follows: 
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At y = 6k(x) 


q>y = ek* (x)cp x (33) 

As in the basic field, so also here the values of the various potential 
coefficients cpr^ on the airfoil surface may he obtained by expansion 
about y = 0: 


cp/^ekUT] = qy^XjO) + ek(x)cp xy ii: (x,0) + . . . 
qy^[x,ek(x)] = cpy^XjO) + ek(x)qp yy lk (x,0) + . . . 


(34) 


Making use of equations (34), substituting equation (6) into equa- 
tion (33)j and ordering in powers of b and 6, the following set of 
boundary conditions results: 

In b: 

<Py 10 (x,0) = -<Py 01 (x,0) = -n(x) (35a) 

In e : 

(Py^XjO) = k* (x)q> x 00 =• k' (x)Mq (35b) 

In b 2 : 

cj^r i:L (x, 0) = -qpyP^x,©) =0 (35c) 

In eb: 

cp y 13 (x,0) = k’ (x) | 9 x 01 (x, 0 ) + q^U/O)] - k(x)|qx yy 01 (x,0) + (p yy 10 (x,0)] 
= k* (x)Jm(x) + cp x 10 (x,02j - k(x)[p 2 m t (x) + <P yy 10 (x,o7| (35d) 


In 


e 


2 . 


c Py 1 ^( x ,°) = k ' (x)cp x 12 (x,0) - k(x)cp yy 12 (x,0) (35e) 


Also in equations (35) the right-hand sides are known after the lower- 
order problem has been solved. 
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Characteristic Condition 

Procedure .- It may be noted that for the computation of the basic 
flow field both velocity components were prescribed, separately, along 
the axis, whereas in the present problem only the ratio of the two veloc- 
ities is fixed by the airfoil condition. Thus, an additional condition 
is required. 

From the fundamental laws of supersonic flow it is evident that no 
perturbation of the basic field can exist upstream of the two character- 
istics passing through the leading edge in a downstream direction. Exist- 
ence theorems for hyperbolic differential equations require that the 
potential be continuous everywhere. Hence, the perturbation potential 
must vanish all along the leading-edge characteristics. 

To formulate this characteristic condition analytically, the equa- 
tion of the characteristics has to be known. Suppose the differential 
equation is given by: 

A(x,y; V^zx + - VV)V + 


C(x,y; (Px^cPyy + F(x,y; cp x ,cpy) =0 (36) 

where A, B, C, and F are linear or nonlinear functions of x, y, 
and the first derivatives cp^ and cpy. Then the equation of the char- 
acteristics is as follows: 

A dy 2 - 2B dx dy + C dx 2 = 0 (37) 

so that 

§ = |( B± 'I® 2 - Ac) (38) 


Thus, the characteristics are seen to depend only on the coefficients 
of the second-order derivatives in equation (36). In general, these 
coefficients are functions not only of x and y but also of the 
desired function cp. Hence, in a rigorous theory the characteristics 
are not known beforehand but only after the function qp is determined. 
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By ass umi ng a power-series solution for cp, this indeterminancy is 
removed; the nonlinear equation, equation ( 36 ), is replaced by a set 
of linear differential equations whose second-order derivatives have 
the same constant coefficients for each equation. In this way, the 
characteristics of each linear equation of the set turn out to he the 
same straight Mach lines, x ± f3y = Constant, determined by the recti- 
linear component of the basic flow field, whereas the true characteris- 
tics are functions of the entire flow field (including the perturbation 
field qp* ) and hence they are curved. This is an unsatisfactory state 
of affairs. 

It is possible to improve the, computation of the characteristics 
somewhat by considering once more nonlinear equation (5), before it is 
split up, and by merely distinguishing between the potential of the basic 
field cpP, which is known, and that of the disturbance field qp 1 , which 
is looked for. If the order of magnitude of the various potential coef- 
ficients qp^ is, for the moment, left out of consideration, curved 
characteristics can be obtained also for linear differential equations, 
as the following illustration may show. 

Substituting the potential series (6) in the form: 



(39) 


equation (5) becomes: 

{1 ♦ i-1-% 2 - (<Py B - V)!] - + <) 2 } (fe B * %.') - 

2(9/ + 0(<Py B + V)(v B * W ) * j^ 1 + " K B + ?* , )3' 

4^ B + V) 2 }(v E + V)-° < to > 

Since cp® has been determined already as a function of x and y, equa- 
tion (4o) may be considered as a differential equation for the unknown cp' 
only and may then be written in the following way: 

ACP^’ + ZBtpxy' + CCPyy' + F ^X , J Cf^' ^y^ = 0 ( ^ 
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Here F is, as before, a given function of x, y, and the first deriv- 
atives <p x ' and cpy* . Each of the coefficients A, B, and C consists 

of one part that depends only on cpB and a second part which is propor- 
tional to qpx' or qpy’ : 



(42) 


Since it is the prime purpose of this analysis to determine the influence 
of the inhomogeneous field upon the pressure distribution around the air- 
foil by means of a linear theory, it is sufficient to approximate the 
equation for the characteristics, by neglecting in the coefficients A, 

B, and C that part which depends on cp' , so that now: 


~ f l( < P B ) ~ gi( x ^y) 


A 

B 

C 


In this way, the characteristics of the total field around the airfoil 
are approximated by the characteristics of the basic inhomogeneous field. 
On the other hand, the lines represented by equation (38) which result 
from such a procedure will be curved. Thus, the equation for the char- 
acteristics becomes: 


jl + - ( c Py B )^] ' ' L T 2 ( c Px B ) 2 J dy2 + 2 K B( Py B ) d y + 

£ + - m 3 - 2 4rH q v B ) s *} = ° (^ 3 ) 


where 
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Analytic formulation .- Substituting equation (44) into equation (43) 
and retaining only terms up to first powers in b, tbe coefficients of 
equation (43) become: 


A = -|j 2 + (r + DMobq^Ol + . . ._ 


2B = -2M o b( Py 01 

C = [l - (7 - l)M o bcp x 01 . . f] 


> (45) 


where cp^l and cpyPl are functions of m(x,y) and n(x,y) given by 
equations ( 24) . 

Inserting equations (45) into equation ( 38 ) and keeping consistently 
only terms up to first powers in b, the equation of the characteristics 
becomes : 


dy _ Mpbcpy 01 ± \)[T - (7 - DMobq^Qg^pg + (7 + DMobq^ 0 ^ 

** P 2 + (7 + l)Mobcp x 01 


(46) 


and after expanding the right-hand side in powers of b: 


Q = ±1 } 

dx P 


1 - ”||E + 1 - - ne!k 01 1 b t ’v 01 • • • } 


(47) 


The upper (or lower) sign is used for the characteristics above (or below) 
the airfoil. 


Equation (47) may be solved by an iteration process, in which the 
first approximation is obviously 
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This leads to the Mach lines: 


= ±x/p 


m 


where the constant of • integration has been set equal to zero so as to 
yield the leading-edge characteristics. 

Substituting equation (48) into the right-hand side of equation (47) 
yields in the second approx im ation: 

I = + x) + (7 ■ 

(49) 

From equations (24) substitute for the characteristics above (or below) 
the airfoil: 


mtSbc) ± + m(0) T 


' tE 

(»y 01 )y=±Vs * lF m(2=t) + ^ T m(0) + ^ 


(50) 


Inserting equations (50) into equation (49) and integrating yield 
for the characteristics: 


y = to(x) + bt^x) . . . 


where 


\lf Q (x) = ±x/p 



The upper (or lower) sign refers to the upper (or lower) surface, 
respectively. 
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Coming back to the characteristic condition indicated in the sec- 
tion "Procedure, " the following requirement has to be fulfilled: For 


y = t 0 (x) + b4r 1 (x) + . . . 

cp’ = bcp-^ + b^cp-*--*- + e<fp- 2 + ebqr^ + 6 2qjl4 = o 


More explicitly: 


bcp 10 f 0 {x) + tnJfpCx) . . 7j + b^cp 11 |x; to( x ) + ^(x^) . . 7] + 

ecp 12 ^; i(r 0 (x) + bijr-^(x) . . 7j + ebqr^jlc; t 0 (x) + bt!(x) . . 7j + 

t 0 ( x ) + ^p(x) - . 7] = 0 (52) 

Expanding cp' about y = ^q(x) in powers of b yields: 

bqr L0 (x, 1 r 0 ) + b 2 f 1 (x) qV 10 (x, ll r 0 ) + cp 11 (x,t 0 )] + e<p 12 (x,i|r 0 ) + 
eb[t 1 (x)qy 12 (x,t 0 ) + cp 13 ( x ^ 0 )] + f V^(x^ 0 ) = 0 (53) 


Splitting equation (53) in powers of b and e furnishes for each 
potential coefficient <p^ one characteristic condition: 

In b: 


c P 10 ( x ^o) = 0 - (54a)' 


In e: 


(p 12 ( x ,\jro) = 0 


(54b) 


In b 2 : 


< P 11 ( x ^o) = -ti( x )9 y 10 (x,t 0 ) 


(54c) 
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In eb: 


In 



, <P 13 (x^ 0 ) = -t!(x)cp y 12 (x^ 0 ) 

qp^fx^o) = 0 


(54d) 

(5te) 


Whenever the relation is inhomogeneous , the unknown potential is expressed 
in terms of known lower-order potentials. 

Each of the boundary-value problems consisting of a differential 
equation, a boundary condition, and a characteristic condition has now 
been formulated and will be solved in the appropriate order. 


COMPUTATION OF FEKTURBATION POTENTIAIS 


In the following discussion (see fig. l) k(x) denotes the ordinate 
of the profile above (or below) the x-axis, which is the direction of the 
relative wind Mq. To distinguish between the upper and lower surfaces, 

k U and kp are used; k(l) is zero or unequal to zero depending on 

whether the airfoil is at zero or at finite angle of attack. If there 
are two signs in front of a term, the upper sign refers to the upper, 
and the lower sign, to the lower surface. 


Computation of 

That part of the disturbance potential which is proportional to b, 
namely, cp^O, obeys the following equations: The differential equation: 


R 2 m 10 m 10 „ r\ 

P 'Pxx ~ Vyy " 0 

the boundary condition: 

cp y 10 (x,0) = -n(x) 
and the characteristic condition: 


(55) 


cp 10 (x,±x/p) = 0 



2k 


NACA TN 2849 


To obtain a solution, it is convenient to introduce the characteristic 
coordinates: 


£ = x + 
n = x - 3y 


Then 


and 


?x = «P g + 'Pn 

= p(q> g - 


(56) 


P 2 <p. 


XX 


- V ■ 


(57) 


As the characteristic condition is different on the upper and lower sur- 
faces, the two surfaces will be treated separately. 

Upper surface .- In characteristic coordinates equations (55) may "be 
written as follows for the upper surface: The differential equation 

becomes 


^ 2 <P ^ 10 = 0 
the boundary condition, 

p ( 9 f 10 - 

and the characteristic condition, 

<p 10 (i,o) = 0 


(58) 


Integration of the differential equation yields: 


<p 10 (6/l) = F(|) + G(tj) . (59) 


where the arbitrary functions F and G are determined by the boundary 
condition and the characteristic condition. 
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From the characteristic condition there follows: 

j 

qP| 10 (^O) = F'(£) =0 for any | 


(60) 


Then the boundary condition requires that 


p[f’(6) - G*(Tij] = p[f*(x) - G’(xj] = -n(x) 


(61) 


Taking into account equation (60), it follows that: 


G’(x) = | n(x) 

p 


G* (tj) = i n(q) 
P 


( 62 ) 


Therefore : 


<P x 10 (x,y) = | n(x - py) 


qpy 10 (x,y) = -n(x - py) 




(63) 


To obtain the pressures on the airfoil surface y = ek(x), equations (24) 
may be expanded: 


<f>x 1 °E J ’ ck U ( x Il = < Px 1 °(^°) + ek u (x)cp xy 10 (x,0) . . 


— n(x) - ek TT (x)n’ (x) + . . . 
P U 


^y 10 E ,ek u(x) = -n ( x ) + e P k u( x ) n '( x ) + . . • 


(64) 


Lower surface.- On the lower surface, the boundary condition is the 


same: 


Kn 10 - V 0 )^ ■ - n(i > 
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but the characteristic condition differs: 


<pl °( X -’"f) = < P 10 ^ 0 ^) = 0 


( 65 ) 


From equation (65) there follows: <p^( 0 ,T|) = G* (q) = 0 , for any t). 

Then_, from the boundary condition there follows: 


3q>£ = PP'U)| =X = -n(x) 


so that 


<Pg(6) = -- n(£) 


( 66 ) 


Using equations (42) 


and 


10/ X 

<Px ( x >y) = 


q>y l0 (x/y) = -n(x + Py) 


n(x + py) 


( 67 ) 


On the airfoil surface: 


<Px 10 jx^ek L (xIl =-i n (x) - ek L (x)n' (x) . . . 
cp y 10 [x, 6 k L (x)] = -n(x) - epk L (x)n’ (x) . . . 


( 68 ) 



HACA TN 28k9 


27 


Computation of qp 1 ^ 

For the potential qpH which is proportional to b 2 , the boundary- 
value problem can be formulated by means of equations (32c), (35c), 
and (54c): The differential equation is 

P 2 ^ 11 - c Pyy 11 = -Wo £|Jr + 1) + (r - i)p^](cp 3OC 01 + qpx X 10 ) (q^ 01 -f-cPx 10 ) + 

2 ( cp y° 1 + < Py 1 °)( < ftty 01 +c Pyx 10 )} + <fyy° 2 - P 2( Pxx° 2 
the boundary condition, >-(69 

cp y i:L (x,0) = 0 . 

and the characteristic condition, 
cp 1;L (x,±x/p) = -ilr 1 (x)qylO( Xj±x /g) 


Upper surface .- Inserting the values for cp 01 , cp 10 , and 

qPyy 02 - P 2 ^^ 02 from equations (24), (64), and (30), for the upper 

surface problem (69), in characteristic coordinates, is as follows: 
The differential equation becomes 



the boundary condition, 

- V 1 )^ - 0 (71) 

and the characteristic condition, 

cp 1IL (5,0) = ♦ 1 (6,0)n(0) 


( 72 ) 
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The rest of the procedure may he indicated here briefly, the details 
being carried out in appendix A. 

Upon integrating equation (70) with regard to £ and T), respectively: 
CD^ 11 = G 1 (m,n) + P-l* U) 

cp^ 11 = G 2 (m,n) + P 2 ‘ (-n) 

where G^ and G£ are known functions of m(g), m(r|), n(£), and 

n(-q) and and P 2 are unkn own integration constants. 

Now, characteristic condition (72), upon differentiating with regard 
to 


(73) 


1,0) = H (!,0)n(0) (74) 

determines P-j_’ (£) by means of equations (73)» 

Then boundary condition ( 71) , after insertion of equations ( 73) > 
fixes P 2 * (t| ) . 

With Pj_ and P 2 - determined, the two velocity components (p x ^(x,0) 
and qpy-*-*-(x,0) that are necessary for the computation of the pressure on 
the airfoil surface can be determined by means of equations (73) and (56). 

Lower surface . - For the lower surface an analogous procedure can be 
carried through; note that the characteristic condition is now: 


^■^-(O,^) = ^r 1 (0,Ti )n(0) 
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The resulting velocity components on either airfoil surface are 
as follows: 

<l V LL ( x ’ 0 ) - 3> M o 2 + 3' ±£ ^ + ± 


n' (x) 


r x n(v_ 

Jn 8 


P _W 0 P 




dt]’ 


+ 2^7 - 1)1^2 + gmUMx). _ 


c Py 11 ( x ^O) = 0 


r (T5) 


Computation of cp^-2 

That part of the disturbance potential which is proportional to e, 
cp 12 , must he expected to equal the expression obtained in the Ackeret 
theory. 

Hence, the velocity components above (or below) the airfoil are: 


<p x 12 - tr fyi) '<* 
9 y 12 = Mo%(l)'^ x + Py) 


(76) 


By expanding about y = 0, the velocities on the upper (or lower) 
airfoil surface are: 


9 x 12& 6k(x y| “*tVl) ,( x) . + * * * 


qyL2|i,dc(xj] “jyfyu'U) T V 6k U(L) (x)k U(L) ' ' (x) . . 


y (77) 
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Confutation of cp-*-3 

The potential (p 1 ^, which is proportional to eh, may be considered 
to originate from some kind of interaction between the inhomogeneity of 
the basic field and the airfoil disturbance. The boundary-value problem 
for qpl3 can be formulated b;y means of equations (32d), (35<i)j and (54d) 
The differential equation is 


- ?>yy 13 " - K 0 {&♦!♦ ' 7 7 1 ) p ll [*“(*,“ + ‘PjOI 10 ) + 
■Pxx 12 ^ 01 * ^ 10 3 - + 'V 10 ) + 

***&* - v°3} 


the boundary condition. 


K78) 


qy>-3(x,0) = %(-£,)’ (x) jm(x) + cp^O^oT] - 


k U(L)( x )[P 2lIl, + <Pyy l0 U,0)] 
and the characteristic condition, 

cp 13 (x,±x/p) = -t 1 (x)cpy 12 (x,±x/p) 


Upper surface .- The treatment of the upper surface is as follows: 

Inserting for cp- 1 - 0 and cp- 1 - 2 expressions (63) and ( 77 ), problem (78) 
written in characteristic coordinates becomes: 

The differential equation: 



( 79 ) 
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The boundary condition: 

3^1 13 - c Pn 13 ) e=11=x = V( x )|m(x) + 2^1] _ p^x) m*(x) + 2^X| (80) 

# 

The characteristic condition: 

cp 13 U,0) = -t 1 (|,0)M o k u *(0) (= 0) (81) 

The rest of the procedure may he indicated here briefly, the details 
being carried out in appendix A. 

Upon integrating equation (79) with regard to | and t|, respectively: 
<Pg 13 = F 1 (m,n,k) + J 1 * (6) (82a) 

and 

Cf^ 13 = F 2 (m,n,k) + J 2 * (q) (82b) 

where Fp and Fp are known functions of k(|), m(E;), n(£), k(q), 

m(Ti), and n(-ri) and Jp and Jg are unknown integration constants. 

Wow characteristic condition (8l) 

q>g 13 (^o) = -ti^^o)^^' (o) 

determines Jp' (S) by means of equations (82). The function Fj 

vanishes for k(0) =0 and k’'(0) = 0, as can be seen from appendix A. 

Also equation (8l) vanishes in that case. Then it follows simply that: 

Jp* (I ) = 0 

This simplification suggests making the assumption, k’ (0) = 0, which is 
an artifice frequently used in linearized theories. The error introduced 
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by this deviation from the true contour is known to yield small inac- 
curacies in the flow field in a very thin strip adjacent to the two 
characteristics passing through the leading edge. Therefore the sim- 
plification is permissible. 

Boundary condition ( 80 ), after insertion of equations (82), fixes 
J 2 * (q). With and J 2 fixed, the two velocity components cp^^O) 

and cpy 1 3( x ,0) that are necessary for the computation of the pressure on 
the airfoil surface can be determined by means of equations ( 82 ) and ( 56 ). 

Lower surface.- For the lower surface an analogous procedure can 
be carried through. Note that now the characteristic condition is 
cft|( 0 ,Ti) = 0 and that different expressions for k(x), cp 10 , and cp 12 

have to be inserted. 


The resulting velocity components on the airfoil surface are: 

<p/3( x , 0 ) = ±|j§ j - 1) + 8 - ^fj k U(L)* ( x ) 

21^2(7 - 3) + 8 + ^|k u(L) (x)[^‘ 


(x) ± 


[»<*> 4 • 

*44 L 


<Py 13 (x,0) = k u ( L )'(x)|m(x) ± _ P 2 k u ( L )(x)jm> (x) ± n ^ j 


(83) 


Computation of cp-*-^ 

O 

That part of the disturbance potential which is proportional to e, 
that is, cpl^, obeys equations (32e), (35e)> and (54e): The differential 

equation is 


P^x^-'Pyy 14 = - M o^+l)cp xx 12 cp x 12 + (7- D^cpyy 12 + 2<p y 12( p X y 12 


the boundary condition. 


(pyl^^O).- qjjj^x^Jk^L)* (x) - q) yy 12 (x, 0 )k u ( L )(x) 
and the characteristic condition, 
qAHx,±x/p) = 0 


H84) 
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A _ _ 

Upper surface. - For the upper surface, after inserting for cp 12 
expressions ( 7T ) > problem (84) , written in characteristic coordinates, 
is as follows: The differential equation becomes 

4(32^^ = - p-[(7 + 1)V (85) 

the boundary condition, 

" cp n l4 )|^ =x = “ T V 2(x) + 3 M 0 k u (x)k u « , (x) ( 86 ) 


and the characteristic condition, 

<p 1 ^(l, 0 ) = 0 


(87) 


Integrating equation ( 85 ) with respect to | and T), respectively: 


9 , 


l4 _ V 

i 4? 




9^ = ■ + 1)%’ (ti)V’^) 0 + Mg'^) 


(88) 


where and Mg have to be determined from the boundary condition 

and the characteristic condition. 


From characteristic condition ( 87 ) there follows: 

9 ^ l 4 U, 0 ) = 0 ( 89 ) 

which, by means of equations ( 88 ), fixes ( 5 ) = 0 if again the assump- 
tion k 1 ( 0 ) =0 is made, as in the computation of cp^-3. 

Then boundary condition ( 86 ), after insertion of equations ( 88 ), 
fixes Mg' (ti ) . 

With Mj_ and Mg fixed the two velocity components on the airfoil 
surface can be determined by means of equations ( 88 ) and ( 56 ). 
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Lower surface . - After analogous considerations for the lower surface 
have been carried out, the velocity components on the two surfaces are: 

<Px l4 ( x >°) = jj(7 + DUq 4 - k (^ 2 - lj] ku(L)’ 2 ( x ) 4 

|^P^u(l) ( x ) k u(L) 1 ' ( x ^^ y (9°) 


c Py 1 ^ x -’°) ~ +LT k U(L)’ 2 ( x ) ± P w o k U(L)^ x ) k U(L)' ’ ^ x ) 


P 


FORCES AMD MOMENTS OR AIRFOIL 
Complete Velocity Field on Airfoil Surface 


To obtain the pressure distribution, the squares of the velocity 
components of the whole field are needed. By adding all of the com- 
ponents (up to the second order in b and e), there result: 

%c 2 " “o' + * 1 * x 01 ) 2 + fe 10 ) 2 + W 0 ] 4 

2 € bcp x 12 (cp x 10 + cp x 01 ) + e 2 (cp x 12 ) 2 + 2M Q [b(cp x 01 + cp x 10 ) + 



= ^(V 2 ) 2 + ^ 2 ( ( Py° 1 + < Py 10 ) 2 + 2ebc ?V 12 (^ 01 4 <Py 10 ) 

Substituting the expressions derived in the preceding sections: 


2 2 2 ? 
<P X = Mq + LA 1 + eAg + b A^ + ebA^ + e A^ 

<J V 2 = e2M o 2k U(L)’ 2 ( x ) 


Y (92) 
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Here: 


r 


n(x) 


A-, = Mq 2m(x) ± 2 

L P 


A 2 = T^-jkudj'U) 


a 3 - 


±a ^] 2 ^^- 3 )t 2 ^]x 


±bM + & (I) ± BlM 


il r sLnLL a ,. 
J J o H 


f> -Wo 

s[Tr - dmo 2 + 2] £wshT] + + i)a£<|i 


M n 3 


k U(L)'( x ) 


y ( 93 ) 


4m 0 2 ( 7 - 1) + 8 - 

VU 


[-<*> * sfcg 


k U(L)( x ) 8 + 2M 0 2 ^ " 3) + ^ P(x) ± £lMl 

L Mo 2 JL p J 


a 5 = - 




Pressure Distribution on Airfoil Surface 

The pressure distribution can be obtained from the velocity distri- 
bution by means of the energy equation (reference 2) which, in dimension- 
less form, is as follows: 



( 94 ) 
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Here p Q , p Q , and Mq are taken for the homogeneous flow field and p, 

p, q, and c are local quantities. By use of the isentropic p - p 
relation, equation (9^0 becomes: 



(95) 


and after substituting p = p Q + Ap and expanding in powers of the 
relative pressure increment — , which is small compared with 1: 


Ap 

1/Al 

’-Y 


— 

6 

II 

1 H 

& 

-lo 

7P 0 

2 \n 

>J • • j 


fo C o 2 2 \po c o 2 / J 


= |j^ 2 - (<P X 2 + <P y 2 ) (96) 


Upon substituting equations (92) for cf^ 2 and qpy 2 , equation ( 96 ) 
may be solved in successive approximations. In the first approximation 
j and second-order quantities in the velocities are neglected. 

Then: 


_^_.. bMo C w ± =fa£l* 

Po c o L p J 





(97) 


Upon substituting from equation (97) 

second-order quantities in the velocities 
second approximation yields: 


for 


(*) 


and keeping 


throughout equation (96), the 



cP-j^ + bP2 + b 2 P2 + ebP^ + s 2 P^ 
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Here 


Mo 


P 1 - ±2 


p 2 - -at|ix) ± 2JiI] 

p 3 - f([r - DHo 4 + 1 —■ - (r + Dm,, 4 2^2>. ± 

2p2 


Ur - 3)^ + 


m' (x) ± 


n l (x) 


r b^ei 

P 


ax’ + 


2(^ 2 - 1 ) 2 «2(x) ± + 1) - K*fe 2 - iJakMiij 

* p ipE r + 1)M ° 4 - H m o 2 - 4] I E (x) 1 ur^uy w + 


f d) ± ! ^]v L )( i '> 


H98) 


p 5 * ^E 7 + ^ - H*fc 2 - ^^(L)' 2 ^) 

2P 


As shown in appendix B, equations (98) for the pressure distribution 
on the airfoil surface can be checked in a limiting case, in which com- 
plete agreement of equations (98) with an expression derived by Ferri 
(reference 2) is obtained. 


Drag Force 

In order to obtain the drag force, the component of the pressure 
which is parallel to the flight direction x has to be integrated over 
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the upper and lower surfaces. Denoting, as in figure 1, by e %(x,) the 

angle subtended by the local tangent of the surface with the flight 
direction, the drag becomes: 


D 


)• 



sin e Sy 


ds 


U 



sin e 6 ds 

J-i Xj 


(99) 


where ds^ and ds L are the line elements of the airfoil surface. 

Keeping only terms up to third orders in s, equation (99) becomes: 


D 


f 


A V 


V (x) 


dx - 



Ap^ek^' (x) dx 


Hence, the drag coefficient is: 




PI 


UO 





dx 


( 100 ) 


For (Ap)y and- (Ap) L expressions ( 98 ) have to be substituted in 
equation (100). 

In the cane of a symmetrical airfoil at zero angle of attack 
ky(x) = -k L (x) = k(x). Then: 



+ bep2 + eb^D^ + e^bDjj. + 
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Here: 



D ^ = _ - 


](7 + 1)M q 4 - 4(^2 - l)| 2 J 1 m(x)k’ 2 (x) dx + 


m' (x)k(x)k l (x) dx 


” 5-^1 


[(7 + UM^ - k^ 2 - l)j J k' 3 (i) dJC 


In the case of a symmetrical airfoil at a finite angle of attack c© 
the decomposition of kg' (x) and k^’ (x), suggested in the section on 
airfoil geometry (appendix C), has to he carried out. In other words: 
Replace V in (Ap)u and k^’ in (Ap) l hy: 


and 


kg' (x) = k(x) - 6 
k^' (x) = -k(x) - Q 


( 102 ) 
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where now k(x) s; k(x) is the customary profile function, for which 
k* (1) = 0. 

Then, upon substituting the modified expressions (98) into equa- 
tion (100), the following generalizations of the coefficients of 
have to be noted: 
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The physical meaning of the various terms in equations (101) 
or (103) may he readily interpreted: (or Dp') is the wave drag 

of the airfoil well-known from the Ackeret theory; (or D^' ) is 

the second-order drag known from the Busemann approximation (refer- 
ence 3 ); Dg represents the horizontal buoyancy acting on the airfoil; 

D 3 may he denoted "second-order buoyancy"; and D]j. is the term charac- 
terizing the interaction between the airfoil field and the given pres- 
sure gradient of the wind tunnel. 

Insofar as bel^ is concerned, it is to be noted that the integral 
is taken over the product of the local volume element ek(x) dx and the 

il . . 

local pressure-gradient coefficient b — m* (x). Since the pressure gra- 

Mq 

dient is varying from point to point, this expression represents, indeed, 
the generalization of the classical result which is thus seen to hold not 
only for incompressible but also for supersonic flow. 

However, the other correction term for the drag, which can be derived 
for incompressible flow (see "Introduction"), does no longer appear in 
the supersonic range: The simple concept of the "apparent mass" is to be 

replaced by more complicated expressions. 


Lift Force 

In order to obtain the lift, the component of the pressure which is 
perpendicular to the flight direction has to be integrated over the upper 
and lower surfaces of the airfoil. Thus: 


L = - J (Ap) n cos s dSy + 


which up to and including second orders 


/(Ap)j 


cos e 


0 L dS L 


in e becomes: 




(io4) 
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The lift coefficient is: 



For (Ap).^ and (Ap) L expressions (98) have to he substituted in 
equation (105). 

In the case of a symmetrical airfoil at zero angle of attack 
kyCx) = -k L (x) = k(x). Then: 


(Cl)^ - “i + b % + sbL ; 


Here: 




(x) 


k(x) dx 


Note that equations (106) vanish for m = n = 0, as they should. 


(105) 


(106) 
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In the case of a symmetrical airfoil at a finite angle of attack F, 
after making the same modifications as in the case of the drag, the fol- 
lowing changes in the coefficients Lj_ have to he made: 


Li* = L,_ 


Lg' = L2 


L 3* " L 3 - ^3 


%V~ 


|Tr + IJMq 1 *' - ^(Mq 2 - l) 2 J m(x) dx 


(107) 


r 1 

/ xm' (x) dx 
J 0 


Furthermore, a coefficient 


appears . 


elV = 


e 


ke_ 

P 


The coefficients L^’ and Lc* represent the lift expressions 
known from the work of Ackeret and Busemann (see reference 3), L-^ 

and 1*2 are vertical forces due mainly to the vertical pressure gradient, 
and is the interaction term. 


Pitching Moment 

The moment about the leading edge is obtained by integrating over 
the moments of the differential lift forces cLL acting at a distance x. 
Therefore the moment coefficient becomes: 



The integral is to be taken over the upper as well as the lower surface. 

Since an upward lift which is to be considered positive produces a 
negative (diving) moment, the contribution of the upper surface is: 



x dx 
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The contribution of the lower surface is: 



In the case of a symmetrical airfoil at zero angle of attack 
kyCx) = -k L (x) = k(x) . Then: 


Ma-o * + b2 “2 + 6bM 3 


Here: 


Mi 


- - l r 1 x » 

Mq J 0 P 


(x) 


dx 


Mo = 


2 2 2 
Mq P 


Hr - 3)% k + <*>* ^T 1 ax ’] 


dx + 


* i) - 2 ( M 0 2 - ^ 


x dx 


m 3 = - 


-i-lHr + l)*, 4 - 4( M o 2 - 41 f^X 1 ^ + 

MqP^ |J- L 0 


n’ (x) 


xk(x) dx 


§ 


(108) 


>(109) 
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In the case of a symmetrical airfoil at a finite angle of attack Q 
after making the same replacements as before: 



Furthermore , two additional coefficients appear: 



and 



* 2 - k < 


E 7 


+ DU, 4 - it(M 0 2 



xk' (x) dx J- 


The coefficients M^' and represent the moments known from 

the Ackeret and Busemann calculations, M 1 and Mg are the moment cor- 
rections due to the pressure gradients of the wind tunnel, and is 

the interaction term. 
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Application of Theory 

As a simple example of the foregoing general considerations, the 
force corrections will he given for a symmetrical airfoil placed at 
zero angle of attack into a tunnel field, for which 


m(x) = n(x) = x 


^i.e., for which the pressure gradient hpUm' (x)-y^- is linear in x^: 


C D - e 


— f fk’ (x7] 2 dx + eh — f k(x) dx + 

0 ^0 L J M o J 0 

e 3 + l)Mo 4 - 4 (Mq 2 - 1 )\J [k'(x 0 3 dx - 
0 0 

eh 2 ^-^2(37 - 17) + 4 (m q 4 + 4 )] J xk(x) dx - 

e 2 h —^-3 [(7 + " ^(Mq 2 " ijj J 1 xk ,2 (x) dx 


«o0- 


- - “ - b2 i^3& 4(57 + 11 - ^ - 2 3 


y ( 111 ) 


Ct — 1) 


eh -i-pR? + 1 )Mq 4 - 4(^2 - l)l r k(x) dx 

Mop 41 - J 0 

° M ■ - b 3V + b£ S?P^ (57 + 1( ■ ^ ' 2 3 

eh — + 1) - 4 (Mo 2 - l)l f xk(x) dx 
M^ 41 - J 0 
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For a doubly symmetric diamond profile of vertex angle e these 
formulas reduce to: 


cd ■ e2 ? 7 si % ■ eh2 Afe 2 ' 37 " 17) + 4 (^ 4 + 4 H 


8(3 


% + x) - 4 ( M ° 2 - 


MqP 




£b ^& 4 < 7 + 11 - ^ 


> ( 112 ) 


C M = -D 


+ 1,2 ^ fc ^ 57 7 7> - 4 ( M <> 2 - £ ll 


3«oP sm^ 


eb + 1} ’ k ^ 2 ' 


8m q P 


Numerical results obtained by means of equations (112) are plotted 
in figures 2 to 9 ‘for representative values of b and e. In figure 2 

the components of Cp which are proportional to e , eb, e b, and eb^ 

are graphed separately, and in figure 3 the total Cp is shown as a 
function of the Mach number Mq. Similarly, the components of Cp and 
Cj^ which are proportional to b, eb, and b^ are plotted separately 
in figures 4 and 6, respectively, while the total Cp and C^ are shown 

as functions of Mq in figures 5 and. 7 , respectively. For the case of 

a finite, small angle of attack eU the behavior of the lift coefficient 
and its components is demonstrated in figures 8 and 9- 
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The relative importance of the various terms may he seen from the 
plots; for example, for the drag the first-order horizontal buoyancy is 
the major correction to be applied and is on the order of 25 percent of 
the true drag. Next in importance is the interaction term which amounts 
to about 10 percent of the wave drag. The second-order buoyancy amounts 
to about 1 percent. 

The same relative magnitude of the correction terms exists in the 
lift coefficient, if the case of an angle of attack is considered as in 
figure 8. 

For the lift and moment at zero angle of attack the term proportional 
to b represents again the largest correction, while the remainder are 
small compared with it. 


New York University 

New York, N. Y. , October 17, 1951 
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APPENDIX A 

COMPUTATION OF PERTURBATION POTENTIAIS 
cp ^ , q>43_, AND cp 1 ^ 

Computation of qpH 


The details of the computation of the potential qpll 
follows : 

Upper surface .- For the upper surface the derivatives 
tions (73) are given by: 


qj. 11 = - «2- ( 7+ i)(i + p2)»(, ) 2hl + 

& 8p 2 ' p 


[(7 + i) + (r - 3)P0 


m*(S) + 


n'(6) 


^ n(y ; 

p 


m(6) + 


1 


PJ P 


y I + iyu) 



are given as 
in equa- 


dq' + 


(Ala) 


S + 


(Alb) 
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Equation (74) in explicit form is as follows: 


1 833 


(7 + 1)(1 + P 2 )[m(0) - si£l + 


[7 + 1 + (7 - 3)P^| 


m 1 


(I) + 


n(j) 


►n(0) (A2) 


The constants of integration are found to "be: 


p i'U) - 

8p 2 1 


(7 + 1)(1 + P 2 )— + 


[7 + 1 + (7 - 3)P^] m 1 ( £) + —Li! 


n(q’ ) 


dq’ 


(A3) 


p 2'(l) = ^5l (7+ 1)(1+ 


n 2 (0) m(q)n(r|) 


P £ 


P 


P P 


>■ + 


n’ (q ) 


{7 + 1 + (7 - 3)P^_, . p 



n(v r j 


m(q* ) + —^3 — dr)' - 




(n) *2lM 


r ^ ^ 

f 0 P 


(A4) 
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Lover surface .- For the lower surface the derivatives corresponding 
to those in equations (73) are: 


^i 11 = (r + 1)(1 + p 2) + m ' T1 + 


E+1+ (7 - 3)pf]-!f^(T,) + 


'-Tf =¥ 9 «' 


y + Pi’d) 


(A5a) 


( P n 11 = + 1) ( 1 + P 2 )m(l)4“ + 


\ 


3 


[r + 1 + (7 - 3)P^ £d) - 


[“' (1) - r *' ^ r| - g-’ 1 at* 


The characteristic condition in explicit form is as follows: 


q^ 11 ^) = ^3 + i)(i + p 2 )j^(0) + ^Tj + 


(A5b) 


[7 + 1 + (7 - 3)pf] jm(ri) - 


nU) 


■n(0) (A6) 
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In this case the integration constants come out to be: 



Computation of cp^-3 

The procedure for computing the potential 'cpl3 is as follows. For 
the upper surface the derivatives in equations (82) are given by: 
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Putting k(0) = k' (0) =0, equation (A9a), together with equa- 
tion (8l ) , yields J-^ 1 (£) =0. On the other hand, substituting equa- 
tions (A9a) and (A9b) into equation (8o) yields: 




(r + i)(i + P 2 ) i 


m’ ( ti) + 




m(q) + 




(n) 


y + 


(7 + 1)(1 + P 2 ) - % 


21 r- 


Mo: 


m(Ti) 




(n ) 


E + 1 + (7 - 3)p5v'(ti) P [m(l’) + dl* + 

7 + 1 + (7 - 3)p 2 + ^2 ^ ^ + n ( A1 °) 


The lower surface can he treated analogously.. 
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Computation of qP-4 

In computing the potential cpl4 the upper surface is treated as 
follows. Putting in the first of equations (88) k(o) = k* (0) =0, this 
equation, together with differentiated equation (87), yields M-^’ (I ) = 0. 

On the other hand, substituting equations (88) into equation (86) yields: 


Mg* (n) = - 



k 2 



(ti) - (r + i)^* (n) 


+ 


k u (n)V' (4) 

V 


(All) 


The lower surface can be treated analogously. 
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APPENDIX B 

CHECK OF EQUATIONS (98) FOR PRESSURE 
DISTRIBUTION IN A LIMITING CASE 

Ferri (reference 2) represents the pressure distribution over an 
arbitrary surface, which is exposed to a uniform stream of Mach num- 
ber M-p as a power series in terms of the local angle of attack tj. 
Thus , 


P - Pi 2 

— = ta^ + agTp + 


p ! c l 


(Bl) 


where a^ and ag are given functions of Mp p is the local pressure, 
and pp the free- stream pressure."*' 

One may note that Ferri* s result corresponds to a limiting case of 
the general nonuniform field, namely, to the case in which the basic 
flow is characterized along the x-axis by an x-component M Q + bm and 

a y-component bn, where both m and n are constant. 

The agreement of the present result with Ferri* s can be seen as 
follows. Expressions (24) for and cpy 0 - 1 - indicate that if m(x) 

and n(x) are constant along the axis, they are also constant in the 
whole xy-plane. The same is no longer true in order b^, because the 

differential equation for cp 02 is inhomogeneous. However, up to 
order b, the basic field may be considered uniform in this limiting 
case. 


The Mach number of this uniform field is: 


M x = — ^(Mq + bm) 2 + (bn) 2 

The angle of the flow direction with the local surf are is 


T| = € 0’ 


+ 


bn 

Mq + bm 


(B2) 


(B3) 


"*"Note that Ferri’ s coefficient is here multiplied with M-, 2 . 
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■bn 

where e0’ and -r- — are the angles subtended by the x-axis with 

+ Dm 

the local surface and with the stream direction, respectively. 

In applying equation (Bl) it should be noted that Ferri's pressure 

P - Pn 

coefficient is taken with respect to M-,, while the coefficient ^2. 

o c 
M o o 

in this report is based upon Mq, which is related to by equa- 

tion (B2). The connection between the two pressure coefficients is 
given by the relation 

P - Pp 

p 

Pl c l 


" Po Po 


c l 2 V p o c o 2 Pl 



(B4) 


where — — , representing the pressure difference between points 

Pl c l 2 

where the Mach numbers Eire Mq and Mp, as well as P 0 |Pp and c o/ c l 
may be obtained by the isentropic energy relation. 

Substituting Ferri 1 s result" for p - p^, with M 2 emd q taken 
from equations (B2) and (B3) and the expression for p Q - p obtained 

from the energy relation, equation (B4) may be expanded in powers 
of b and yields: 

P - Po B '% 2 / n \ 

* 2 ^ - $ - 

* pr((7 + 1 > M o lt - k ( M 0 2 - l)] (m ± |) (E5) 

p 

which includes terms of order e, b, sb , and 6 but not terms of 
order b^. 
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On the other hand, equation (B5) can he obtained from equation (98), 
if m(x) = m, n(x) = n, and k' (x) = -S’ are substituted and the 
b 2 term is disregarded. Thus, in the limiting case considered, Ferri* s 
and the present results are identical. 



58 


NACA TIT 28^9 


APPENDIX C 
AIRFOIL GEOMETRY 


In order to show that the considerations in the text pertain to 
airfoils at both zero and finite angles of attack, the geometry of the 
airfoil will he considered in detail. 

If the chord of the airfoil is at zero angle of attack, the local 
slope is given by 


— = sk* (x) £ 60 (Cl) 

dx 


where y and k(x) are positive on the upper, and negative on the 
lower surface (see fig. 1(a)). If the airfoil is symmetrical, then: 



ek* (x) ~ e Qjj 



-ek* (x) s: e0 L 


(C2) 


Figure l(b) may illustrate hcrw these relations have to be modified 
in case the chord of the airfoil is at a finite angle e0 relative to 
the flight direction x. 


Two coordinate systems may be used, one whose x-axis coincides with 
the flight direction and another one whose x-axis coincides with the 
chord. Then the slope of the profile above the flight direction x is: 


dy 

dx 


€k* (x) 


= tan e(0 - 6) 


6(0 - 0 ) 


(C3) 
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where tan 60 is the slope of the profile above the chord. If the 
profile function above (or below) the chord, y = ek(x), is introduced, 
then 


60 Z tan e0 = k* (x) 


and equation (C3) can be rewritten: 



(c4) 


Note that y > 0 on the upper surface and y < 0 on the lower surface. 
Then, at P in figure l(b), 60 is_positive because y increases_with 

increasing x. Fixing the sign of 0 analogously to that of 0, 0 in 
the sketch is positive and larger than 0, so that in equation (C3) ^ 
comes out negative, as it should according to the figure. 

Now, the barred and unbarred coordinate systems are related (neg- 
lecting second orders in e) as 'follows: 

x ~ x - y0 
y £ y + x 0 

Then 

ek' (x) = ek' |x(l + - 0 k(x)e^j 

~ ek 1 (x) 


Therefore, equation (c4) may be written: 


= gk' (x) = g[k ’ (x) - 0~1 
dx . ■— — 1 


(C5) 
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Figure .1.- Airfoil geometry. 
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Figure 4.- Components of lift coefficient plotted against Mach 

number Mq for diamond profile. € = 0.05; h - 0.05. 
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Figure 5-- Total lift coefficient 
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Cj, as function of Mach number Mq 
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Figure 8.- Component of lift coefficient plotted against Mach 

number Mq foi* dismcnd. profile si^t qttiqI i sjiglo of 

attack €0. 60 =. 0.05; e = 0.05; "b = 0.05* 
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